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Abstract 

Within the Extended Nambu Jona-Lasinio model, we calculate the 0{p^) countert¬ 
erms entering the low-energy expansion of the 77 —> and the rj —>■ 7 r °77 am¬ 
plitudes in Chiral Perturbation Theory. For 77 —our results are compatible 
with both the experimental data and the two-loop calculation using meson reso¬ 
nance saturation. For the p decay we find F (7 —>■ Tr^yy) = 0.58 ± 0.3 eV which is in 
agreement with experiment within one standard deviation. We also give predictions 
for the neutral pion polarizabilities and compare them with the results obtained 
from resonance saturation. 
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1 Introduction 


Over the last few years the transitions 77 —> and 7 —>■ 7 r °77 have been stnd- 

ied in the framework of Chiral Pertnrbation Theory (yPT). These two processes 
have in common the fact that their amplitnde starts at 0(p‘^) and that there is no 
contribntion from tree diagrams at this lowest order. 

The first transition has recently been analysed in [|l| at the two-loop level. The 
0 (p®) connterterms have been evalnated assnming resonance satnration and it is 
shown that the predicted cross-section fits well the existing data [Q]. It has been 
noted by several anthors p -0 that the contribntion of these connterterms to the 
cross-section is very small in the threshold region. On the other hand, when not 
only high-precision data from DA<hNE become available for energies np to 600 MeV, 
bnt also the theoretical predictions are improved by nnitarization procednre (this 
reqnires calcnlating the fnll 77 — 7 r’'' 7 r“ amplitnde to two-loops), it may become 
possible to extract one linear combination of the 0 (p®) conpling constants. 

Concerning rj a particnlar featnre of this decay is that the 0{p'^) one-loop 

contribntion to the width is more than two orders of magnitnde smaller than the 
experimental measnrement. In this snppression has been explained in physical 
terms and a partial analysis of higher loop contribntions has been carried ont. It is 
expected that the contribntion of the 0{p^) connterterms shonid acconnt for a large 
part of the amplitnde. The analysis of gives a reasonable estimate of the order 
of magnitnde for the decay width. However the latter is still a factor two smaller 
than the experimental valne [§| r{p vr^yy) = 0.85 ± 0.18 eV. New data may 
become available in the foreseeable fntnre at SATURNE in Saclay, where a proposal 
for measnring this decay width has been recently approved. 

The S-matrix elements for both transitions depend on the same set of 0{p^) 
conpling constants. Thronghont the following we restrict onrselves to the large Nc- 
limit, where only three low-energy conplings enter these transitions. The theoretical 
analysis is not yet refined enongh and the present experiments are not snfficiently 
precise, to extract the valne of these conpling constants from the data. Therefore it 
is important to make predictions for these conplings (this is particnlarly relevant in 
the case of the p vr'^yy decay width which has a strong dependence on the conn¬ 
terterms). One possibility is to estimate their valne nsing the resonance satnration 
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method, which provides successful predictions to 0{p'^) [P, |^. However it is not 
known whether the 0{p^) couplings are actually saturated by resonance exchange. 
In addition to the estimate based on resonance saturation, one can hnd in the lit¬ 
erature a calculation of two of the 0{p^) coupling constants using the chiral quark 
model [0. 

The purpose of our paper is to compute the three 0{p^) coupling constants of 
the large iVc-limit in the framework of the Extended Nambu Jona-Lasinio (ENJL) 
model. Within this model the low-energy effective action of QCD has been derived 


n| to O(p^), as well as some of the coupling constants governing the low-energy be¬ 


haviour of the lowest spin-1 and spin-0 resonances ||11|, The model has only three 
parameters (which is an advantage compared to the resonance saturation method 
where the number of parameters increases when new transitions and higher chiral 
order are considered) and the agreement with experiment is good. 

This paper is organized as follows. In section 2 we present the structure of the 
0{p^) chiral Lagrangian in the large iVc-hmit. Section 3 contains the expressions for 
the amplitudes of the transitions and the pion polarizabilities. Section 4 is devoted 
to the ENJL predictions for the 0(p®) coupling constants. Finally in section 5 we 
discuss the numerical results and compare our predictions with those of the resonance 
saturation method and with the existing data. 


2 0{p^) chiral Lagrangian 


Let us collect as usual the octet of Goldstone bosons in the unitary unimodular 
matrix U 


U = exp —i\/2 


fo 


( 1 ) 


where /q ~ /tt = 93.2 MeV and (A are Gell-Mann’s SU{3) matrices with trXaXb = 
2Sab) 
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The axial-vector field matrix is defined as follows: 


- ^[df, - (4) 

where = U and are external 3x3 left and right field matrices. We also 

define 


/(7) = ei"re^±e^i"re, (5) 

where are the external field-strength tensors 

^ ^ ^ ( 6 ) 

Since we are only concerned with /(_,_) we set in what follows = f^'^. The 
specification to the electromagnetic field reads 

= \e\QF>^‘^ (7) 

where Q = diag(2/3, —1/3, —1/3) is the quark charge matrix. Finally we shall need 

( 8 ) 

with 

X = 25oAf, (9) 

where A4 = diag(m„, nid, mg) and Bq is related to the vacuum expectation value 




( 10 ) 



The structure of the strong chiral Lagrangian up to has been studied in 


two neutral pseudoscalar mesons and two photons. By restricting ourselves to the 
operators whose coupling constants are leading in the large-iVc limit we can write 

= 4 + d2tT (e„e/.0/“'’) + 4tr {x*ur ^)}. (n) 

J 77 

We have considered here that all quantities were communting since we are interested 
only in neutral particles. Our conventions and notation are chosen in such a way 
that the expansion of this Lagrangian in terms of pseudoscalar helds coincides with 
the analogous expansion of the Lagrangian (4.28) in |]^ which was introduced for 
the SU{2)l X SU{2)ji case. 

In the following section we display how the physical amplitudes are expressed in 
terms of these three coupling constants di, d 2 and d^. 


13 . It is straightforward to write the relevant 0(p°) chiral Lagrangian involving 


3 Amplitudes in xPT 

The Lorentz covariant and gauge invariant amplitudes for 7 (q'i) 7 (q' 2 ) —^ 7 r°(pi) 7 r°(p 2 ) 
and r]{pi) n^{p 2 hiqi)'l{q 2 ) read 

A = e^e^(qi)e^(q2)V^^, (12) 

where 

= A(s, f, + B(s, f, u)T 2 f,^, 

rr 

J-ifiu c\9iJ.v qiuq2^j,i 

2 (13) 

= 2 sA^A^ - {t - - 2(t - u){qi^/A^ - g 2 /.A^), 

= (pi -P2)^l■ 

3.1 77 ^ TT^TT® 

For the 77 —>■ amplitude the Lagrangian o generates the 0{p^) counterterms 
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10 


= -^d: + 


(14) 


We refer to Ref. []^ for the calculation in SU{2)l x SU{2)r of the additional con¬ 
tributions coming from chiral loops indicated here by the ellipses. The contribution 
of one kaon loop in the 0{p^) amplitude has been computed in [jl^. For s far below 
the KK threshold it is numerically small, with respect to the contribution of one 

I^ , ||. The calculation of the two-loop 


pion loop, owing to an extra factor 


48mi 


amplitude in SU{3)l x SU{3)r has not been carried out. 

As in [|^ one connects the constants a^, and 6^ parametrizing the renormalized 
amplitudes 


with the constants df 


aN _ 
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3.2 T] 

For the p decay we have [0 

A^/2 


(16) 


^6 = 3^"4 [16(dgm^ - dlmi) + {d[ + 8d^)s] + 

_ — + ■ ■ ■ . 

3v/2/4 


(17) 


In these expressions we have included the contribution of ds which has not been 
taken into account in ||^. However the size of this contribution is small because 
<< m^. The mixing rj — rj' has been treated within yPT in [^, |^. As a result 
the physical r/ is a superposition 
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rj = cos Oris — sin Orji , 


(18) 


with sin 0 ~ 

For the expressions of the O(p^) one-loop contribution to this process we refer 
to Ref. [0. As explained there, this contribution is very small because pion loops 
violate G-parity and kaon loops are suppressed by a factor l/24m|'. The two-loop 
calculation has never been performed but the same argument can be advocated 
1]^, to claim that loop contributions to the 0(p®) amplitude are suppressed. In 
addition, consistency of the yPT expansion implies that they are smaller than the 
0{p‘^) one-loop contribution. 

On the other hand at 0{p^) a new type of one-loop contribution appears by 
taking two anomalous vertices of 0{p‘^). As shown in [^] the latter has an order 
of magnitude comparable with the 0{p‘^) one-loop contribution. This fact does not 
break the perturbative yPT expansion because the higher order corrections to this 
0{p^) one-loop term will be small with respect to it. 

Because all these loop contributions are small with respect to the measured decay 
width, one should expect that the 0{p^) counterterms account for a large part of 
the full amplitude. 


3.3 Neutral pion polarizabilities 


The polarizabilities caracterize the electric and magnetic properties of a composite 
system. They appear as parameters in the low-energy expansion of the Compton 


amplitudes at threshold [17 


jiCompton _ 2 


ei ■ ^2 


(— - a^a;ia; 2 ) - (gi x ei) ■ (gs x es*) H- 


, (19) 


where gf = (cuj, ^). Following the Condon-Short ley phase convention we dehne 


(a + f3) = 8 Q;m 7 rlim lim B, 


a 


( 20 ) 


[a — 13) = — lim lim (A -|- 8m^B), 
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The one-loop amplitude calculated in has been used in Eqs. in order to 

hnd the neutral pion polarizabilities to 0{p^), as discussed in Refs. |IB| , |^. Taking 
into account the full 0{p^) result one gets 


(« - Pf = 

{a + (3)^ = 


+ 44 - 44 ) + (« - 


80 


—am^. 


—d\ + (a + ^)f_ioop 


where the two-loop contributions can be obtained from Table 3 in [p] 


( 21 ) 


( 22 ) 


(“ - /3)2-loop - -0-31. 

(« + /3)f-loop=^0'17. 

The latter numerical values have been derived for the SU ( 2)4 x SU{2)r case; we are 
neglecting the additional two-loop contributions for the SU{3)l x SU{3)r case due 
to kaon loops. Notice that the kaon loop contribution to the one-loop amplitude, 
calculated in |T^, vanishes for s —> 0. Hence, at the one-loop level, there is no kaon 
loop contribution to the pion polarizabilities, according to (^). 


4 ENJL model prediction for the coupling con¬ 
stants 

The idea of the ENJL model is to approximate large-iVc QCD at the chiral symmetry 
breaking scale by an effective four-fermion theory (generated by integration over 
quark and gluon helds above A^). The main assumption of this model is that higher 
dimension fermionic operators are irrelevant for long-distances. The three parame¬ 
ters are the scale A^ and the two coupling constants Gs and Gy of the four-fermion 
operators, respectively the scalar-pseudoscalar and the vector-axial ones. Alterna¬ 
tively one can trade these parameters for three other ones, i.e. the constituent quark 
mass, Mq, the coupling of the constituent quarks to the pseudoscalar current, gA, 
and the ratio Mq/A"^. 

*The values of the polarizabilities are in units of 10“"^fm^ in what follows. 


7 











Integration over quarksQ below (see |H] for the description of the method) 
yields the low-energy effective Lagrangians involving pseudoscalar mesons and the 
lowest spin-0 and spin-1 resonance. Within this context twenty-two low-energy 
constants have been derived in and thirty-six others in 


12 . Whenever there 


are experimental data, the ENJL prediction is in good agreement with them. 

A calculation of the coupling constants di and d 2 in the context of the chiral quark 
model which is nothing but the mean held approximation of the ENJL model, has 
been carried out in Ref. (in this article however the operator modulated by d^ 
in Eq. O has not been taken into account). In this mean held approximation one 
neglects the huctuations of the resonances so that the calculation involves just one 
loop of a constituent quark of mass Mq. 

In the full ENJL model (see Ref. [^) the quark-loop contribution is modihed 
with respect to the mean-held approximation by the presence of the mixing axial- 
pseudoscalar parametrized by the constant qa- In addition one gets contribution 
from integrating out the resonance helds. For any coupling our notation is (in 
what follows we shall drop the index ‘r’ of the coupling constants) di = di + 

R 

where di is the contribution of the quark loop and stands for the sum of the 

R 

contributions of the resonances {R = S, V, A, T). 

Calculations have been performed with the Seeley-DeWitt expansion [^, ^ . For 


details of the procedure the reader can consult Ref. 0 . A brief summary is given 
in appendix. 

For all the numerical applications we shall use for the input parameters the values 
obtained in [O, Fit 1, i.e. Mn = 265 MeV, qa = 0.61 and x = = 0.052. 


4.1 Quark-loop contribution 

The results read 


tHere we disregard the gluonic fluctuations below (see for a discussion of this issue). 



















(23) 


where we have used the shortened notation r„ = r(n, Mq/A^) for the incomplete 
Gamma functions dehned as follows: 



(24) 


^ JMyKl Z 

For applications we use Fi = r 2 = 1. When one sets gA = ^ our results for di ^2 
coincide with those of the quark-loop calculation of [Q. 

4.2 Spin-1 and spin-2 resonances 

Here we shall follow the notations and dehnitions of 0 . 0 . 0 - The vector (re¬ 
spectively axial-vector) resonance held will be denoted by (respectively A^j. We 
also dehne 

Rfiu = - d^Rf,, for R = V, A with d;, = + F^ , 

where F^ is given by the expression 




(25) 


Let us hrst consider the vector case. The coupling we are interested in is the one 
which modulates the operator e^^^'^tr(V)i{^y,/po-}) called hy- The 0(p®) countert¬ 
erms receive contributions from the exchange of a;°, and 0 mesons. The last one 
will not be considered here since its contribution is suppressed by the large 0 mass 


(the 0 resonance contribution to the 77 —>■ amplitude has been included in 


Refs. P, m). Guided by the nonet assumption, we shall not use for the contribution 
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of the meson the experimental data on the decay —>■ which have a large 

error, but instead those on the decay Then from r(p+ —>■ vr+y) and 

r(a;° —>• Try) |§] one can extract 


h h 

^ ~ 4.7-10“®MeV^ and ^ ~ 4.9 ■ 10“® MeV 




M^o 


(26) 


and take, according to the nonet assumption, 

“pO “p+ 

Within the ENJL model the constant hy has been computed in |^. Here we 
just recall the result 




(27) 


where fv modulates the operator tT{f^uV^'') and has the following expression |]n 



Nc 2 
167r2 3 


Using the relation found in the ENJL model 
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one obtains 


fvMv 


U 

VI - 9a 


(28) 


(29) 


hy 


iVc V2 

167r2 8/^ 


(1 + 9a)\/^ 


- gA- 


Numerically one gets, for = 0.61, 


(30) 


ft 


V 


Ml 


= 3.6 ■ IQ-^MeV 


7-1 


(31) 


V 


For this constant the agreement of the ENJL prediction with experiment (^) is not 
as impressive as for the other low-energy constants. Moreover this quantity has to 
be squared in the contribution to the amplitude, so that hnally the ENJL model 
prediction of is a factor two smaller than the one based on resonance saturation. 
However one has to keep in mind that the coupling constants predicted in the ENJL 
model are parameters of the Green functions evaluated at zero momenta. This fact 
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has to be taken into account, when comparing with processes, such as the decay of 
a vector resonance, whose typical energy scale is large. In [Q it was argued that 
including intermediate resonance exchanges and chiral loops one obtains an improved 
prediction in good agreement with the phenomenology of the vector resonance decay. 

Whenever we deal with a processes at small energy, such as Compton scattering 
at threshold for the determination of the pion polarizabilities, we favour the value 
of hv predicted in the ENJL model. On the other hand, for the p decay the energy 
scale is quite close to the p mass and it seems appropriate to use for hy the value 
extracted from experiments. Finally for the 77 —>■ transition we shall display 

two different cross-sections corresponding to the values of hy obtained from ENJL 
and phenomenology, respectively. We take the difference between the cross-sections 
as a contribution to the uncertainty on our result. 

The h\ axial-vector resonance is not included in the ENJL model. We shall use 
the experimental data in order to incorporate it as it is done in Q]. 

The tensor resonance is not described either by the ENJL model. Such a reso¬ 
nance would appear only if higher derivative four-quark operators were taken into 
account. In this case again we shall use the experimental data. Unfortunately they 
allow only for a determination of the absolute value of the tensor contribution to the 
transitions considered here. Hence this contribution will appear as an uncertainty 
in our results. 


4.3 Scalar resonance 


The scalar-sector analysis requires the knowledge of the coupling constants C 5 , 
C™ and C 5 which modulate respectively the operators tr(S'.^^.^'^), tr(S'x’'') and 
. The coupling constants as well as the scalar mass Ms 


have been computed within the ENJL model in |^. Their expressions reads 




= 


Nc Mq 
Ibvr^ Xs 


‘^9 A (ro — hi). 

(32) 

p (r-i — 2ro), 

(33) 


with the rescaling factor 
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and the mass 


N 2 

A| = —^-(Sho -2ri) 

^ 167r2 3^ ^ 


(34) 


M| = 


Nc 

Ihvr^ A| ° 


(35) 


Numerical evaluations for these constants can be found in as well as the com¬ 
parison with experiment and with the resonance saturation approach. 

We have computed the coupling constant Cg which governs the decay S' —> 77 
with the result 



Nc 2 hi 

167r2 3 MqXs ' 


(36) 


Numerically we hnd Cg ~ 0.19 GeV“^. This value is about a factor 2 larger than the 
estimate given in [Q (N.B. with an uncertainty of 100%). However it corresponds 
rather well to the estimate obtained in the analysis of ||^ , where higher order terms 
in the chiral expansion are taken into account, in order to match the QCD high 
energy behaviour. 

The contributions of the scalar particle to the constants di (or 01^2 and b) are 
given in |I| in terms of the mass Ms and the absolute values of the constants Cg, 
C^i C's (in contrast with the resonance saturation calculation, here we know the 
sign of the different contributions). Within the ENJL model we hnd 
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o' 


chcifi 

( Nc 4 \ 

1 

8 M| 


' 48 To 

ClC^fl 

( Nc 4 \ 

1 

00 

1 


' 96 To 


gl (U - 


p(r-i - 


r.) 


2r„). 


(37) 


(38) 


(39) 


We hnd af ~ 0.14 and af ~ 4.7 whereas the resonance saturation approach used in 
Ref. 0 gives af ~ ±0.8 and af ~ ±1.3. 
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A few comments are in order. First we have a definite prediction for the signs, 
even though it is very dependent on the input parameters in the case of af which 
is proportional to the difference of df and df, two quantities of comparable order 
of magnitude. Secondly and as a consequence of the latter remark we conhrm the 
smallness of af. Finally we get for af (which is proportional to df) a sensibly higher 
value. This is related to the fact that in the ENJL model the constant Cg computed 
in Eq. (|36D appears to be higher than the estimate of []^ (but, as we said, this is 
in agreement with the recent analysis of ||^) and the scalar mass rather low (see 
|]). At the present time it is not completely clear if the scalar contribution to 
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the coupling constants predicted by the ENJL model is a better estimate than the 
analogous contribution calculated in [|^ . Hence we keep the latter as a lower bound 
estimate and obtain in this way the uncertainty on the scalar contribution to the 
coupling constants. 

In the next section we discuss a possible way of improving the situation in the 
future which is related to the coupling constant ds. 


4.4 On the constant 

An interesting feature of the scalar sector is that within the ENJL model many 
cancellations occur between the quark-loop and the scalar resonance contributions. 
This is the case at 0{p‘^) for the constants L 5 , Lg and H 2 (see [^). In our calcula¬ 
tion, one hnds that the combination af J-di is actually independent of the parameter 
p and of the quadratically divergent incomplete gamma function F_i. For the full 
dg one obtains a quite simple expression 


dg — df -|- dg — 


' Nc fl \ 1 
Ibvr^ Mq j 96 




(40) 


Let us also notice that we get the following relation^ 


Cl 


dg 




(41) 


With the set of parameters we have already used one finds from Eq. (^^ dg ~ 
1.5 X 10“® whereas the value predicted in n using scalar resonance saturation is 


^which is valid to all orders in the asW-expansion (see Ref. |^T)) 
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ds — 0.38 X 10“®. The difference between these values conies both from df (which is 
already bigger in our case than the value of []^) and from wich is not present in 
the resonance saturation approach. One can compare both values with the estimate 
(is = (0.94 ± 0.47) X 10“^ obtained in using a sum rule for the vector-vector 
two-point function. This estimate is affected by a large uncertainty wich makes it 
compatible with both our value and the value in [||] and so does not provide a very 
satisfactory test. 

An important improvement would be to extend this sum rule to two-loop order, 
as it is required by the fact that ds modulates an 0(p®) operator. Hence it is a bit 
premature to try to improve the agreement of our result with the sum-rule, which 
could be achieved by changing slightly the value of the parameters Mq and qa and 
redoing the fits of Ref. including ds] this would of course have an effect on all the 
0(p®) constants. We postpone such considerations until a full two-loop calculation 
of the sum-rule is carried out. 


5 Results 

5.1 77 ^ TT^TT*^ 

In the Table we display the full ENJL results for the constants ai, 02 and b to be 
compared with Table 2 in [^. 

Table 1: Final results 



Q.L. 

(1--) 

(1+-) 

(0++) 

(2++) 

Total 



ENJL 

Data 


ENJL 

Data 



ai 

-12.3 

-20.3 

-36.6 

0 

0.14 

0.8 

74.1 

-28 > ai > -53 

02 

6.1 

7.6 

13.7 

-1.3 

4.7 

1.3 

±1.0 

13 < 02 < 24 

b 

1.0 

1.3 

2.3 

0.7 

0 

0 

±0.5 

2 < 6 < 4 


Inspection of the Table shows that we have three sources of uncertainties. The hrst 
two are related to the different estimates of the vector and scalar coupling constants 
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obtained in section 4 by using ENJL and experimental data, respectively]^. In the last 
column we have added these uncertainties in such a way as to have an upper bound 
and a lower bound for the cross-section. This is because the uncertainties on the 
vector (see section 4.2) and the scalar resonance (see section 4.3) contribution have 
different origins. The last uncertainty comes from the sign of the tensor contribution 
which is undetermined. We included also this contribution in the total (last column 
of the Table). 

For the low-energy constants entering the helicity amplitudes we get 

—8 > = ai -|- 86 > —17 , 

8 < h, = 02 - 26 < 15 , (42) 

2 < = 6 < 4 . 

Our results are compatible with those quoted in |l|: 

h+ = -14 ± 5 , 

h, = 7 ± 3 , (43) 

= 3 ± 1 . 

In the Figure we display the data from Ref. for the cross-section (t(s; |cos 6 *| < 
Z = 0.8) as a function of the center-of-mass (c.m.) energy E = ^/s. The full 
(respectively, dashed) curve in the Figure corresponds to the values in the l.h.s. 
(respectively, the r.h.s.) of the inequalities (|42|) . For comparison we show in the 
Figure, as a dotted curve, also the cross-section obtained in for the central values 
of (|43|) (see the full lines in Figs. 5,11 of Ref. [Q). We can also compare our results 
with Fig. 9 of Ref. |l|], where the uncertainties in the values of (^3]) are taken into 
account (in Fig. 9 of Ref. |I|, however, the contribution of the integrals Aa,b coming 
from the two-loop box and the acnode diagrams in Fig. 4 of Ref. is neglected, 
whereas in our Figure the full contribution is retained). The upper curve in Fig. 9 
§We assume the central value of the experimental data, ignoring the errors on the measured 
parameters (which typically have a 10% size). In addition, we disregard the changes in the ENJL 
prediction due to variations of the three parameters of the model. In our opinion, assuming 
the quite large band of error on the cross-section associated to the ENJL prediction and the 
experimental estimate for the resonance contributions, is pessimistic enough. 
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of Ref. [m disagrees with our result by more than 20% for c.m. energies E >540 
MeV (at E = 540 MeV the ^ contribution to the cross-section is very small, i.e. 
1.5%, see Figs. 5,9 of [|l|). 

5.2 T] 7r®77 

As we said in section 3.2 the 0{p^) counterterms contribution is expected to account 
for a large part of the amplitude because the loop contributions are suppressed. In 
Ref. [0 it is shown that the resonance saturation approximation is not sufficient to 
explain the experimental decay rate. Indeed if one takes into account only vector 
resonance saturation one gets 


r(r 7 —>• 7r°77) ~ 0.18 eV, (44) 

to be compared with the experimental value r (7 —>■ 7r°77) ~ 0.85 ± 0.18 eV. Let 
us review the calculation of |^. By keeping the momenta dependence in the vector 
propagator the authors of 0 made an ‘all-order’ estimate, i.e. at 0{p^) and higher, 
of the counterterms contribution with the result 


T(r] —> 7r°77) ~ 0.31 eV. (45) 

Then if in addition one takes into account chiral loops (the 0{p‘^) ones as well as 
the 0(p®) ‘doubly-anomalous’ one-loop diagrams) one reaches the value 


r (?7 —>■ 7r*^77) = 0.42 ± 0.20 eV, 


(46) 


where the error includes scalar and tensor contributions (whose sign was not known) 
as well as a 30% error coming from other contributions such as one-loop diagrams 
involving the 0{p^) Lagrangian (pd]). 

Let us now show step by step the ENJL results. To the estimate (^4|) which 
incoporates only vector resonance exchange, we hrst add quark loop contribution 
and hnd 


T{p —>■ 7r°77) ~ 0.36 eV. 


(47) 
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Already we are not far from the central value of Eq. (|46|) . In addition if we include 
the scalar contribution as predicted by ENJL we have 

r(77 —> 7r°77) ~ 0.5 eV. (48) 

We shall use as lower bound the prediction for the scalar exchange given in |]l[ . Then 
we have 


r(? 7 —>• 7r°77) = 0.45 ± 0.05 eV. (49) 

Finally taking into account the axial resonance and, as in |^, the contributions from 
chiral loops and, in addition, an uncertainty from the contribution of the tensor 
resonance, we get 


r (7 ^ 7r°77) = 0.58 ± 0.12 eV. 


(50) 


Of course we can play the game of making an ‘all-order’ estimate of the counterterm 
given by vector resonance exchange as it is done in |^. Starting from the upper 
bound in the last equation, we obtain a decay rate of 0.86 eV. However this represents 
only a partial resummation of the chiral corrections. First there are 0(p®) one-loop 
corrections involving the 0{p^) Lagrangian O which could be as important as the 
countertems contributions (even though they are next-to-leading in 1/W)- Secondly, 
in the large W-hmit, other chiral corrections {i.e. of 0{p^) and higher) coming from 
the counterterms occur within the ENJL model. Indeed the quark loop contribution 
considered here corresponds to the computation of a four-point function at zero 
momenta. The full ENJL calculation at non-zero momenta (as it has been done in 


2 ^ for the two-point and some of the three-point functions) would bring additional 


chiral corrections. We shall include all chiral corrections of this type in an additional 
30% uncertainty. Our prediction is then 


r( 7 —i> 7r°77) = 0.58 ± 0.3 eV. (51) 

As we said above, the higher order corrections included in the ‘all-order’ estimate 
of the vector resonance contribution enhance the decay width and therefore go in the 
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right direction, in order to match the experimental result r(r 7 —7r°77) = 0.85±0.18 
eV. It would be interesting to carry out the 0(p®) loop analysis (this can be done, 
given that most of the loop diagrams are suppressed (see section 3.2)), as well as to 
calculate in the ENJL model the 0(p®) counterterms. In this way one could obtain 
an accurate description of the 0{p^) chiral corrections. 


5.3 Polarizabilities 


The pion polarizabilities have been estimated through dispersion sum rules [17 


{a-/3)^ = -10±4 , 

{a + f3)^ = 1.04 ±0.07 . (52) 


At order 0(p®) resonance saturation gives 

{a-(5)^ = -1.01-0.31-0.58 ±0.20 =-1.90 ±0.20 , 

(a + /3)^ = 0.00 ±0.17 ±1.00 ±0.30 = 1.17 ±0.30 , (53) 


where the three contributions that add up to the hnal results on the r.h.s. are the 
one-loop result, the two-loop contribution, and the one from the 0(p®) counterterms 
assuming resonance saturation, respectively. There is a large contribution from 
the 0(p®) low-energy constants due mainly to the a;-exchange in the resonance 
saturation approach For the estimate of the uncertainties in (|53D the reader is 
invited to consult Ref. 0]. 

The contribution due to the 0(p®) counterterms to the pion polarizabilities 


(« - P)e,. 

{a ± /?)^4 


arrin ^ 


am-„ 

(47r/^)^ 


8h_ , 


(54) 


is obtained from the ENJL model prediction for the low-energy constants displayed 
in Eqs. (|42| ) 


-0.33 >{a- > -0.69 , 

0.78 < (a ± /3)^i < 1.44 . 


(55) 
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Adding the one-loop and two-loop contributions - identical to those in - we get 
our results with the ENJL model 


-1.65 > > -2.01 , 

0.95 < {a + f3)^ < 1.61 , 


(56) 


which are compatible with Eqs. ([53|). Both the prediction in (|56|) and the result in 
m are compatible with the forward sum rule {a + j3)^ = 1.04 ± 0.07 in (^). 

As we mentioned in section 4.2, we consider the ENJL value of the vector 
resonance coupling as favoured. This choice corresponds to the range of values 
— 1.65 > (a — > —1.68 and 0.95 < (a - 1 - < 1.28. Here the residual uncer¬ 

tainty for {a — j3)^ reflects an incomplete knowledge (see section 4.3) in estimating 
the coupling constants for the scalar resonance, whereas the undetermined sign of 
the tensor contribution to the low-energy constants in the Table yields the remain¬ 
ing uncertainty on (a -|- {3)^. The sizeable difference between the value of (a — {3)^ 
obtained using the ENJL prediction for the vector resonance coupling and the result 
(|5^ is mainly due to the fact that the ENJL prediction of the vector contribution 
to h_|_ is a factor two smaller than the one based on resonance saturation and the 
experimental data (see section 4.2). The size of the quark-loop contribution, which 
is included in the ENJL prediction, is not large enough, to bring the ENJL value 
for to the one obtained from resonance saturation. 

The construction of unitarized S'-wave amplitudes for 77 —tttt which contain 


(a — (3)^’^ as adjustable parameters has been carried out in Ref. |^. In this case. 


only {a — j3)^'^ can be determined from the data [Q, with the result 


[a- pf = 4.8 ± 1.0 , 

{a-(3)^ = -1.1 ±1.7 . 


(57) 


The value (^7]) for [a — j3)^ is consistent with the two-loop result for the neutral 
pion, both for the calculation based on resonance saturation and for the ENJL 
prediction, whereas the corresponding two-loop calculation for charged pions is not 
available and so it cannot be compared with the value (pTl) for (a — (3)^. 

In the validity of the errors quoted in a recent estimate of (a ± (3)^’^ by 


Kaloshin and collaborators is questioned. Here the polarizabilities appear as 
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adjustable parameters in the unitarized D-wave amplitudes, hence the values of 
(a + can be determined from the data with the result 


+ = 0.22 ±0.06 [28] , 

(a + /3)^ = 1.00 ±0.05 [| . 


(58) 


The authors of arguing on the partial wave analysis of the data that shows 
large uncertainties even at the /2(1270) mass, conclude that the errors quoted in 
(|5^ for {a ± f3)^ are unbelievably small. This result is compatible with the ENJL 
prediction (^), as well as with the value in (^). 


6 Conclusion 

Within the ENJL model we have computed the 0{p^) coupling constants entering 
the xPT expansion for the amplitudes of 77 —>■ vr^vr'’ and 7 —>■ 7 r° 77 . In addition to 
the contribution of the resonance exchange, one has to add terms coming from the 
constituent quark loop. Both contributions have comparable orders of magnitude. 


This situation at 0{p^) differs from the one at which is described in [IT . 

Our major result concerns the decay 7 —>■ where the contribution of the 

counterterms dominates the amplitude. The values of the couplings calculated in 
the ENJL model yield a prediction for the decay width which is compatible with 
the data. Our central value compares much more favourably with the experimental 
result than the value calculated in within resonance saturation. 

Concerning the transition 77 —our results are compatible with those of 
the resonance saturation approach worked out by [|l[], as well as with the data 0 ]. 
For the neutral pion polarizabilities, we hnd results consistent with the estimates 
given in earlier references [Q, 

To improve our knowledge of the coupling constants dj’s and be able to test in a 
more accurate way the ENJL prediction, we need both experimental and theoretical 
progress. 

High-precision data from DA<hNE on 77 — 7 r° 7 r° may allow to extract the value 
of hs, since the cross-section in Fig. 9 of [^ shows a sizeable dependence on this 
low-energy constant for energies near 600 MeV. It will be necessary to carry out a 


unitarization of the two-loop result, using a procedure analogous to [31] and match- 
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ing the dispersive calculation for the all-order amplitude with the two-loop ampli¬ 
tudes for 77 — 7 r° 7 r° and 77 — , once the latter will have been calculated. 
The consideration of this improved, unitarized amplitude will justify the inclusion 
in the analysis of experimental data up to 600 MeV. 

For the other source of information, i.e. the decay 7 —it will be very 
interesting to investigate the 0{p^) analysis, as we said in section 5.2. The crucial 
test may be the future SATURNE experiments which may allow to determine di, 
d 2 as well as the size of the 0{p^) chiral corrections. 

Finally the determination of the constant d^ from the sum-rule has to be improved 
by carrying out the computation to the two-loop order. 
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Appendix 


Let us define 


, 

= U - |75(s/.& + w'-i) - , 

Let us also introduce the following expression: 

M = -MQ-A-i(S-75A) , 

where 

A = . 

Here S, HA+) and denote the helds of the resonances, respectively scalar, 

vector and axial-vector. 

We recall briefly the Seeley-DeWitt expansion of the operator 

~ Mq = — + E , 

where Ve is the Euclidean Dirac operator 

Ee = IfN fj. + M 


and 


E = 2MqS + 


(here we do not need the complete expression which can be found in O)- 
The effective action reads then 

1 dr ^ \ ^ Nc ^ Tin —2, x) 

2 J 1 /AI T E xj 2167r2^o 
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in terms of the incomplete Gamma functions defined in (M). We give the Seeley- 
Dewitt coefficients up to a total derivative and a circular permutation. Only the 
terms we are interested in are displayed here. These terms read 


H3 — —- — -ES/fj^E + —'V^R^p'VyRup ( , 


-^^4 — ^ + — [ERpyERpy + AE‘^ RpyRpy] 


We have also made use of the following identity: 


^^{{dpfupf] = ti{2[{dpfpyf + ifpufupfpp + ipfypfppiu] 
1 
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